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SUMMARY 


This paper describes work currently in progress to update an existing 
transonic circulation control airfoil analysis method. Existing methods 
suffer from two deficiencies: -the inability to predict the shock structure 
of the underexpanded supersonic jets; and the insensitivity of the calcula- 
tion to small changes in the Coanda surface geometry. A method developed 
for the analysis of jet exhaust plumes in supersonic flow is being modified 
for the case of the underexpanded wall jet. In the subsonic case, the same 
wall jet model has been modified to include the calculation of the normal 
pressure gradient. This model is currently being coupled with the transonic 
circulation control airfoil analysis. 
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theoretical APPROACH 


• POTENTIAL FLOW 

- JAMESON'S FULL POTENTIAL METHOD 


• VISCOUS FLOW 

- WALL JET REGION 

FINITE-DIFFERENCE METHOD USING AN EDDY VISCOSITY MODEL 
FOR CLOSURE 

- AIRFOIL UPPER SURFACE AHEAD OF SLOT AND ALL OF LOWER SURFACE 

INTEGRAL METHODS 

LAMINAR -- COHEN AND RESHOTKO 
TURBULENT -- GREEN ET AL. 


Figure 2. TRACON--Component Analysis Procedures. 
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INPUT AIRFOIL GEOMETRY, 
ANGLE-0 F-ATTACK, 
BLOWING MOMENTUM 
COEFFICIENT AND 
REYNOLDS NUMBER 


COMPUTE POTENTIAL 
FLOW SOLUTION 


CALCULATE VISCOUS 
LAYER PROPERTIES 
AND SEPARATION PRESSURE 

P sepu AND P sepl 


DETERMINE SOURCE 
DISTRIBUTION REPRESENTING 
DISPLACEMENT AND 
ENTRAINMENT EFFECTS 
OF VISCOUS LAYERS 


ESTABLISH NEW ESTIMATE 
OF CIRCULATION 
AROUND AIRFOIL 


COMPUTE POTENTIAL FLOW 
ABOUT CONFIGURATION 
INCLUDING VISCOUS EFFECTS 


■—YES- 
1 ' 


PRINT RESULTS 


(Psepu " P sep| ) ** 0 
AND 

\ C L n “ C Ln_i > 


STOP 


Figure 3.. Calculation Procedure. 
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Having obtained the points of separation and corresponding 
static pressures for the upper and lower surfaces, a new value of 
circulation is estimated on the basis of separation pressures and 
current value of lift. A new potential flow solution can then be 
computed using this new estimate of circulation with the viscous 
effect, i.e., velocity component normal to the surface taken into 
account. 

Convergence is checked at this stage. The calculation con- 
tinues for another cycle unless P S epu and P S epl are in close 
agreement and the variation of the lift coefficient between 
successive iterations is in the range of convergence. 

As shown in Figures 4 and 5, the basic analysis performed 
quite well. A relatively high Mach number case (0.6) at an angle 
of attack of -10° produces a strong shock wave as shown in Figure 
4. The analysis accurately predicts the shock location and 
overall pressure distribution. Similarly in Figure 5 for a low 
Mach number but with blowing, the analysis again performs well in 
comparison with experiment. The strengths and weaknesses of 
TRACON are summarized in Figure 6. 


3 . 0 CURRENT DEVELOPMENTS 

The two basic objectives of the current work are as follows. 
(1) Formulation of the methodology to describe the inviscid flow 
of an underexpanded wall jet into a co-flowing stream. The 
method must be able to describe the supersonic flow that results 
from a nozzle operating choked into a stream having a free-stream 
Mach number in the range 0.3 <L M 0.8. As part of this objec- 
tive, procedures for coupling the wall jet calculation with the 
overall circulation control calculation method, TRACON, were 
considered. (2) Formulation of methodology to determine the 
effect of the Coanda jet on surface pressures. In this approach, 
the direct influence of the normal momentum equation is required. 
As with the first objective, procedures must be developed for 
coupling the normal pressure gradient calculation with solutions 
for the streamwise momentum and with the external potential flow 
(TRACON) . 

A literature survey indicated that although several re- 
searchers have investigated the flow field of supersonic jet 
plumes, including the underexpanded jet, only one group, from 
SAIC/Princeton led by Dash and co-workers (5), has considered the 
underexpanded wall jet. The shock structure of an underexpanded 
wall jet is physically very similar to that of a jet plume; 
because of this. Dash and Wolf (6) have recently modified one of 
their codes for application to curved wall jets. For the viscous 
model the finite difference, two-layer turbulence model of Dash 
and Beddini (7) has the required theoretical basis for applica- 
tion to circulation control wall jets. This analysis procedure, 
called SPLITWJET (hereinafter termed WJET) , is completely des- 
cribed in (7) and (8). WJET solves a set of curvilinear, higher- 
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o a DATA; NCCR 1610-8054S 

a = - 10 * 

M = 0.6 

C 4 = 0 

TRrtCON 



(CIRCLE = UPPER SURFACE) 
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SUMMARY 


* ZERO BLOWING — SHOCK LOCATIONS PREDICTED QUITE 
ACCURATELY FOR ROUNDED TRA I L I NG-EDGE AIRFOILS 


* LOW BLOWING - PREDICTED PRESSURES AND C L - C u 
IN GOOD AGREEMENT WITH EXPERIMENT 

* HIGH BLOWING — PREDICTED BEHAVIOR GOOD AT LOW 
MACH NUMBERS BUT AT HIGH MACH NUMBERS A MODEL OF 
THE JET SHOCK STRUCTURE IS NEEDED 

* BECAUSE OF EMPIRICAL MODEL TO ACCOUNT FOR COANDA 
SUCTION, THE EFFECT OF SMALL CHANGES IN GEOMETRY 
CANNOT BE ACCOUNTED FOR IN A CONSISTENT MANNER 

* DISPLACEMENT THICKNESS COUPLING BETWEEN VISCID AND 
INVISCID FLOWS MAY BE INADEQUATE AT HIGHER MACH 
NUMBERS 


Figure 6. TRACON — Summary of Program Features. 
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order boundary layer equations employing a mapped surface— normal 
grid network. A key ingredient of WJET is the inclusion of a 
hybrid, two-layer turbulence model which couples a damped Van 
Driest inner mixing layer solution to a two-equation, ke outer 
solution, matching at the grid point corresponding to y+ ~ 50 (y+ 
= nr s ipi/Ujj,), where n is the distance from the surface, t s is 
the wall shear stress, p is the density and is the laminar 
viscosity. The turbulence model contains curvature corrections 
based on the work of Launder et al. (9) and Hah and 
Lakshminarayana (10) which are required to account for the sig- 
nificantly enhanced wall jet mixing rates associated with strong 
convex curvature. The combined system of mean flow and turbu- 
lence model equations are integrated using an efficient upwind 
fully implicit difference algorithm, and the jet growth is con- 
trolled by an ordinary differential equation keyed to the edge 
vorticity. 

The inclusion of normal pressure gradient terms due to 
curvature in a boundary layer approach has been addressed by 
Mahgoub and Bradshaw (11). For the wall jet application the 
following approach has been developed. 

1 * Matching between the potential flow and the wall jet occurs 
at the wall jet edge. 

2. The pressure variation across the wall jet is solved using 
the viscous normal momentum equation. 

Since the inclusion of the normal momentum equation renders 
the wall jet equations elliptic, direct spatial marching of the 
complete set of wall jet equations is ill posed. The remedy 
involves splitting the pressure field such that: 

1. the wall jet equations ( not including the normal momentum 
equation) are solved in a parabolic fashion with the stream- 
wise pressure gradient determined from a globally imposed 
pressure field, P(s,n); and 

2. a revised pressure field, P(s,n), is determined in the 
course of the parabolic wall jet solution using the con- 
tinuity and normal momentum equations; the revised pressure 
field is used in local thermodynamic relations and the 
normal velocity distribution arrived at from the coupled 
continuity/normal momentum equation solution at each step is 
used in the convective terms of the parabolic wall jet 
equations. 

This procedure is summarized in Figure 7. 

In the subsonic case the strategy for coupling WJET with 
TRACON (Figure 8) is as follows: 
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SUBSONIC WALL JET CALCULATION 



* PARABOLIC SOLUTION WITH d/ds [P(£ ,77)] FROM POTENTIAL 
FLOW SOLUTION 


* CROSS-FLOW SOLUTION BEGINNING AT OUTER EDGE OF WALL JET 
WITH BOUNDARY CONDITIONS FROM TRACON SOLUTION GIVES NEW 

p 

SECOND ITERATION WITH d/ds FROM UPDATED 

PRESSURE FIELD 

Figure 7. Analysis Procedure for Subsonic Wall Jet. 
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OVERLAP COUPLING FOR A SUBSONIC WALL JET 



* SOLVE POTENTIAL FLOW OVER TRACON MATCHING SURFACE WITH 
<^ n (s) PRESCRIBED 

* SOLVE WALL JET USING PRESSURE SPLIT PNS APPROACH 

* DETERMINE 5*(s) FOR OVERLAP REGION AND GENERATE ^ n (s) 
FOR NEXT POTENTIAL FLOW SOLUTION 

Figure 8. Coupling Procedure for Subsonic Wall Jet. 
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!• the TRACON potential flow calculation is performed utilizing 
a trailing-edge geometry which smoothly blends the slot lip 
to the airfoil trailing edge (this surface is defined as the 
TRACON matching surface in Figure 8; and 

2. the WJET wall jet calculation is performed over the actual 
airfoil surface (Figure 8) . 

To interactively couple the two solutions, the WJET calculation 
requires an inviscid flow map of the TRACON calculation to obtain 
outer edge conditions and streamwise pressure gradients, while 
the TRACON calculation requires the stipulation of the source 
strength distribution, <f>n(s), along the TRACON matching surface. 
The definition of the source strength distribution in a flow with 
large normal pressure variations requires careful derivation of 
displacement thickness relations. The usual simplifying approxi- 
mations employed for standard boundary layer problems cannot be 
employed. To date, the wall jet analysis has been coupled with 
TRACON and preliminary calculations are being performed to check 
out the program logic. 

In the supersonic (under expanded wall jet) case a similar 
strategy has been developed. In recent years considerable re- 
search has been carried out on the analysis of supersonic plumes. 
Dash and co-workers at SAIC/Pr inceton (12) through (14), and 
Wilmoth (15) of NASA Langley Research Center have developed shock 
capturing technology which accurately models the multi-cell em- 
bedded shock structure of under expanded jets. Steger (16) and 
Diewert (17) of NASA Ames research Center, Birch of Boeing (18), 
Perry and Forrester (19), Shang of AFWAL (20) and Cline of Los 
Alamos (21) have all in recent years considered the analysis of 
supersonic free jets. These analyses have ranged from two- 
dimensional parabolized Navier-Stokes solvers to full three- 
dimensional time-dependent codes. The shock wave structure of an 
underexpanded wall jet is very similar to that of a jet plume, 
with one exception: the wall jet does not have the extended 

inviscid core region that characterises underexpanded free jets, 
and hence, the waves in wall jets would be propagating in a fully 
turbulent environment several slot heights downstream of the slot 
exit, see Figure 9. To properly treat this viscous problem, a 
parabolized Navier-Stokes formulation is required which can be 
directly coupled with a potential flow solution. 

The basic structure for a supersonic wall jet solver already 
exists in WJET. In order to avoid having two separate codes to 
analyze subsonic and supersonic wall jets, it is more practical 
to develop a supersonic pressure solver for WJET. Thus, WJET 
will be used in its present form to solve the entire wall jet 
with: 

1. pressures in supersonic regions determined by a character- 
istic based pressure solver; 
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SUPERSONIC WALL JET CALCULATION 



* SUPERSONIC RE6I0N — PARABOLIC SOLUTION WITH d?/ds FROM 
VISCOUS CHARACTERISTIC METHOD (WAVE SOLVER) 

* SUBSONIC REGION -- PARABOLIC SOLUTION WITH <?P/<9s FROM 
POTENTIAL FLOW SOLUTION 

* SONIC LINE MATCHING AS IN FREE JET APPROACH OF DASH 

* CROSS-FLOW SOLUTION IN SUBSONIC REGION WITH BOUNDARY CONDITIONS: 

V(77 = 0) - u < 77 = 0) TAN0 CHARCi 
PI.TJ- 1) = P e 

Figure 9. Analysis Procedure for Underexpanded Supersonic Wall Jet. 
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2. pressures in subsonic regions split with the streamwise 
gradient imposed from TRACON and the normal variation deter- 
mined by solving the continuity and normal momentum equa- 
tions; and 

3. coupling procedures utilized to interface the pressures at 
the near wall jet and mixing layer sonic lines. 

The same overlap procedure described previously will be 
utilized to couple the wall jet and TRACON potential flow solu- 
tions. However, the matching surface cannot be arbitrarily 
blended from the slot lip to the airfoil surface as if it is for 
the subsonic case, since it must remain above the jet mixing 
layer sonic line (Figure 10). The extension of the overlap 
approach to accommodate this modification is conceptually 
straightforward, but the details require significant considera- 
tion. 


4 . 0 REFERENCES 


1. Dvorak, F.A. and Choi, D.H., "Analysis of Circulation-Con- 
trolled Airfoils in Transonic Flow", JL*. Aircraft . Vol. 20, 
No. 4, April 1983. 


2 . 


Jameson, A., "Numerical Computation of Transonic Flows with 
Shock Waves", International Union q . f Theoretical and Applied 
Mechanics . Springer-Verlag, New York, Inc., September 1975, 
pp. 384-414. 


3. Brune, G.W. and Manke, J.W., "An Improved Version of the NASA 
Lockheed Multi-Element Airfoil Analysis Computer Program", 
NASA CR-154323 . March 1978, pp. 69-87. 

4. Green, J.E., Weeks, D.J. and Brooman, J.W.F., "Prediction of 
Turbulent Boundary Layers and Wakes in Compresible Flow by a 
Lag-Entrainment Method", Royal Aircraft Establishment TR- 
72231, December 1972. 


5. Dash, S.M. and Sinha, N., "Noniterative Cross-Flow Integra- 
tion Procedure for the Pressure-Split Analysis of Two-Di- 
mensional, Subsonic Mixing Layer Problems", To be Published 
in the AIAA Journal. 


6. Dash, S.M., Beddini, R.A., Wolf, D.E. and Sinha, N., "Vis- 
cous/Inviscid Analysis of Curved Sub- or Supersonic Wall 
Jets", AIAA Paper No. 83-1679, Danvers, MA, July 1983. 

7. Dash, S.M. and Beddini, R.A., "Viscous/Inviscid Analysis of 

Curved Wall Jets: Part 2, Viscous Pressure-Split Model 

(SPLITWJET) ", Science Applications, Inc., Princeton, N.J., 
TR-7 , November 1982. 

8. Dash, S.M. and Sinha, N., "Pressure-Split Extensions of 
SPLITWJET Model for Wall Jet/Potential Flow Coupling", 


179 


OVERLAP COUPLING FOR A SUPERSONIC WALL JET 



* MATCHING SURFACE CONFIGURED TO RESIDE ABOVE SONIC LINE 

* POTENTIAL FLOW OVERLAP INCLUDES ONLY SUBSONIC REGION 
OF JET 

* P(s) — SUPERSONIC (77 — 0) - WAKE SOLVER 

— SUBSONIC (77 — 0) — PRESSURE SPLIT 

* 5*(s) DETERMINED FOR SUBSONIC PORTION OF JET ~ INCLUDES 
INFLUENCE OF JET UNDEREXPANSION 


Figure 10. Coupling Procedure for Underexpanded Supersonic Wall Jet. 
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